(N 



The z/MSM, Dark Matter and Baryon Asymmetry of the Universe 

Takehiko AsakgQ and Mikhail Shaposhniko\ 1 
Institut de Theorie des Phenomenes Physiques, 
Ecole Polytechnique Federale de Lausanne, CH-1015 Lausanne, Switzerland 

(Dated: May 2, 2005) 

^ - Abstract 

■ We show that the extension of the standard model by three right-handed neutrinos with masses 



smaller than the electroweak scale (the z/MSM) can explain simultaneously dark matter and baryon 
asymmetry of the universe and be consistent with the experiments on neutrino oscillations. Several 
constraints on the parameters of the i^MSM are derived. 
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Introduction. — The canonical Minimal Standard Model (MSM) despite being ex- 
tremely successful in particle physics, cannot accommodate experimental data on neutrino 
oscillations simply because neutrinos are exactly massless in the MSM and thus do not 
oscillate. In addition, the MSM does not contain any particle physics candidate for cold 
dark matter and cannot explain the baryon asymmetry of the universe. 

The simplest renormalisable extension of the standard model, consistent with neutrino 
experiments, contains A/" right-handed SU(2)xU(l) singlet neutrinos Nj {I = 1, . . . ,A/') with 
the most general gauge-invariant interactions described by the Lagrangian: 

SC = Niid^YNi - F^i L^Nj^ - ^ NfNj + h.c. , (1) 

where $ and La (a = e, /i, r) are respectively the Higgs and lepton doublets, and both Dirac 
(Md = and Majorana (Mj) masses for neutrinos are introduced. We have taken a 

basis in which the mass matrices of charged leptons and right-handed neutrinos are real and 
diagonal, and F is a matrix with elements Fd. 

In addition to quite a large number of dimensionless Yukawa couplings, this model con- 
tains Af dimensionful parameters - the Majorana masses of right-handed fermions. The 
neutrino oscillation experiments cannot fix these new scales, as the masses and mixing an- 
gles of active neutrinos contain only specific combinations of Md and Mj, coming from the 
diagonalization of the complete mass matrix. 

In this paper we propose to choose these unknown mass parameters to be of the order 
of the electroweak scale or below. In other words, we will assume that the model defined in 
is a true low energy theory up to the Planck (or, say, grand-unified) scale. Moreover, 
we fix A/" to be 3, keeping the number of right-handed neutrinos to be equal to the number 
of fermionic generations. The specific model in this parameter range will be called below 
the "z/MSM", underlying its minimal character and the fact that no new energy scale is 
introduced. 

The aim of the present work is to demonstrate that the z/MSM with a particular choice 
of parameters, consistent with the data on neutrino masses and mixing, can explain simul- 
taneously the dark matter and baryon asymmetry of the universe. 
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General properties of the uMSM. — The z/MSM contains 18 new parameters in comparison 
with the MSM. Three of them are the Majorana masses, while another 15 are hidden in the 
Yukawa matrix F^j and can be chosen as 3 diagonal Yukawa couplings, 6 mixing angels and 
6 CP-violating phases in the following way: 



F = KLfdKR^ 



(2) 



where 



/d = diag(/i,/2,/3), KL = KLPa, Ka^ = Ka^Pp 



(3) 



The diagonal matrices for Majorana phases are 



P„ = diag(e^"i , e^"^ , 1) , Pp = diag(e^'^i , e'^^ , 1) , 



(4) 



and the KM-like mixing matrix Kl is given by 



Kj 
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1 

CL23 

-SL23 CL23 
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-SLize'^^ CL13 



^ ( CL12 SL12 ^ 



-SL12 CL12 

1 



(5) 



where c^ij = cos^Ouj) and SLij = sm^Ouj)- The expression for Kji is written in full analogy 
with the replacement of L to R. We fix the indices in such a way that "1" corresponds to 
the lightest right-handed neutrino and "3" to the heaviest. The new light neutral fermions 
can naturally be called the "dark" {Nd), for the dark matter candidate, "clear" (Nc) and 
"bright" (Nb) sterile neutrinos. 

As the Majorana masses are assumed to be of the order of the electroweak scale or below, 
the model can only be consistent with the neutrino experiments if the Yukawa couplings are 
very small, ~ 0{mi,Mj/v'^), where rrii, are the masses of active neutrinos and v = 174 GeV 
is the VEV of the Higgs field. At the same time, in the interesting parameter range discussed 
below, the Majorana masses are much larger than the Dirac masses, so that the see-saw 
formula Jsj for the masses of the active neutrinos is applicable. The active neutrino mixing 
matrix is coming then from the diagonalization of the see-saw mass. 
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The vMSM and dark matter. — For small /j the lifetime of a right-handed neutrino may 
exceed the age of the universe; in this case the z/MSM provides a particle physics candidate 
for the warm dark matter P]-!^- The allowed mass range of the corresponding sterile 
neutrino is severely restricted as 

2 keV < < 5 keV , (6) 

where the lower bound comes from the cosmic microwave background and the matter power 
spectrum inferred from Lyman-a forest data [l^, and the upper bound is given by the 
radiative decays of sterile neutrinos in dark matter halos limited by X-ray observations . 

In the recent paper fl3] we have shown that the minimal number of sterile neutrinos, 
which can explain the dark matter in the universe, is A/" = 3. In this case only one sterile 
neutrino can be the dark matter. We identify it with the lightest one A^^i, i.e. A^^i is the 
"dark" sterile neutrino. In addition, the observed mass density of the dark matter leads to 
the following constraint on the parameters of the model: 

{m\,Md)ii ^ ml (7) 
where mo = (9(0.1) eV. In terms of Yukawa couplings and mixing angles it reads: 

/i ^\.i2 '^ms + fl '^ms '^m2 + fl -^ms — 3-3 ■ 10 (8) 

In other words, at least one of the couplings fi must be of order 6 ■ 10^^'^. We choose it to be 
fl. Though being much smaller than the Yukawa constants in the charged lepton sector, this 
value does not contradict to anything and is stable against loop corrections in the z/MSM. 

The dark matter constraint allows us to fix the absolute values of the active neutrino 
masses as follows: the mass of the lightest active neutrino should lie in the range mi < m^™ = 
C(10-5) eV, while other masses are mg = [4.8l|]:^] " lO'^ eV and ma = [9.05^1^;^] ■ lO^^eV 
([4.7^05] ■ 10~^eV) in the normal (inverted) mass hierarchy jl^ . 

There are further constraints on Yukawa couplings coming from Big Bang Nucleosynthesis 
(BBN) and from the absolute values of neutrino masses. Since the second and third active 
neutrinos are much heavier than the first, one must have /2,3 ^ /i. This means that these 
neutrinos equilibrate before BBN and spoil its predictions, unless they decay before BBN. 
This leads to the constraint [3| ^2,3 > 1 GeV, and correspondingly, to /| > 3 • 10^^^, /| > 
2 ■ 10~^^, where the specific numbers were derived assuming small mixing angles Onj and 



normal hierarchy. Together with eq. (jH)), this tells that the mixing angles and Ojn^ 
must be small, with 

^.n < 3.3 . 10- (^) aud < 1.4 ■ 10- (^) . (9) 

For the case of inverted hierarchy the conclusion on the smallness of sru and sris is the 
same. 

We would now like to see whether this model can account for the baryon asymmetry of 
the universe in a specified parameter range. 

The uMSM and baryon asymmetry. — In considering the problem of the baryon asym- 
metry we take the most conservative point of view and assume the validity of the standard 
Big Bang theory well above the electroweak scale. Since the right-handed neutrino Yukawa 
coupling constants are small, A^'s are out of thermal equilibrium at high temperatures and 
may come into thermal equilibrium at smaller temperatures. In any case, the lightest sterile 
neutrino, playing the role of dark matter, never equilibrates. Because of this we will assume 
that initial concentrations of right-handed neutrinos are equal to zero. 

To the best of our knowledge there has been only one study of baryon production in 
the model (P) with small Majorana masses^. Namely, in ref. Akhmedov, Rubakov and 
Smirnov (ARS in what follows) proposed an interesting idea that the baryon asymmetry 
can be generated through CP-violating sterile neutrino oscillations. For small Majorana 
masses the total lepton number of the system, defined as the lepton number of active neu- 
trinos plus the total helicity of sterile neutrinos, is conserved and equal to zero during the 
universe's evolution. However, because of oscillations the lepton number of active neutri- 
nos becomes different from zero and gets transferred to the baryon number due to rapid 
sphaleron transitions Roughly speaking, the resulting baryon asymmetry is equal to 



the lepton asymmetry at the sphaleron freeze-out [IQ], with the exact relation in |17|. 

According to the ARS computation, the produced lepton asymmetry is entirely expressed 
through the Majorana masses Mj of sterile neutrinos, Yukawa couplings fi and the mixing 
matrix and can easily be of the required order of magnitude. However, the dark matter 
and BBN constraints put severe limitations on the phase space of the i/MSM and we found 



^ Note that the Majorana neutrinos with masses much larger than the electroweak scale give rise to thermal 
leptogenesis Q- 
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that no choice of the parameters, consistent with ()8l9p . can lead to the observed baryon 
asymmetry of the universe, if the ARS equations are used. This can be understood in the 
bllowing way. The CP-violating effects must be proportional to the Jarlskog determinant 
isl ] related to the matrix Kji, J = €^1^1^011120^233 mss ji^s R23 sin5/j. With the constraints on 
the mixing angles discussed above, the value of J is at most 10"^''. Other factors, such as 
the total number of degrees of freedom, make the asymmetry predicted by ARS formula well 
below the observed value in the region of parameter space we are interested in. However, 
we have reached a different conclusion. Below we will reconsider the neutrino oscillations in 
the early universe and identify the similarities and crucial differences with ARS. 

In general terms, the evolution of possible lepton asymmetries can be found with the use 
of kinetic equations for a complete neutrino density matrix ^ - [^|. In our case this is a 
12 X 12 matrix with components describing the mixing of active neutrinos and anti-neutrinos 
as well as the sterile neutrinos of different helicity states: 

/ Ot T Ott Ot at Or at \ 



P 



Pll Pll Pln Pln 
Pll Pll Pln Pln 
Pnl Pnl Pnn Pnn 
\ Pnl Pnl Pnn Pnn / 



(10) 



where different entries in p are the 3x3 matrices describing the neutrino states in different 
sectors (by N we denote the negative helicity state of the sterile neutrinos). This density 



matrix satisfies the kinetic equation 



2l[ which can be written in the form 



(11) 



where H = k{t) + H^ + Hint is the Hermitian effective Hamiltonian incorporating the medium 
effects on neutrino propagation [2^, k{t) ~ T is the neutrino momentum, T and are the 
destruction and production rates correspondingly. Following ARS, we will use the Boltzmann 
statistics and replace the last term in (fTT|) by zF^. 

If the system is in thermal equilibrium, the equilibrium density matrix, given by p'^'^ = 
exp {—H/T) must satisfy eq. (fTTj) . which gives the relation F^ = ^{T,p'^'^} ^. 



^ For numerical estimates p'^'^ can be safely replaced by exp {—k/T), as was done in ARS analysis, since the 
corrections to p'^'^ coming from H'^ and Hint are small. 
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In the leading approximation (all Yukawa couplings fi are neglected) the Hamiltonian 
if ° is diagonal and can be written as 

= diag(ii°^, Hl^, HU, H'^^) (12) 

with, in turn 2311 . 



" kit) [ 32 



(13) 



where the first term comes from the electroweak gauge correction to the active neutrino 
propagator, and the second from the Higgs correction, h^, /i^ and are the charged lepton 
Yukawa couplings. For the sterile neutrinos we have 

^^;v = H'^N = ^diag(Ml^ Ml Ml) . (14) 

Because of the structure of the consideration of oscillating neutrinos can be simplified 
significantly. To remove the trivial time dependence of the density matrix, make a change 



P 



U{t)pU\t) , (15) 



where U = exp (— i J^dt'H^{t')), which brings us to the evolution equation in the interac- 
tion picture. Now, as we will see later (and in accordance with ARS) the most important 
region for leptogenesis is Tl ~ (AM^Mq)^, where AM^ is the typical value of the Majorana 
mass squares and Mq ~ 7 ■ 10^'' GeV appears in the time-temperature relation, t = 
At this time most of the off-diagonal elements of density matrix p undergo very rapid os- 
cillations because of U{t) and decouple from the system. These elements are listed below: 
Pln, Pln, Pln, Pln, Pnl, Pnl, Pnl, and p^i; they can be safely put to zero. In ad- 
dition, all non-diagonal parts of pll and pn can be neglected as well. In physics terms, 
the oscillations between L and N are strongly suppressed since these particles have very 
different effective masses. The same is true for La —>■ Lp, a ^ (3 oscillations. On the other 
hand, the non-diagonal parts oi p^N and p^vTV must be kept, as the corresponding expo- 
nentials in U{t) are of order 1. Another four non-diagonal entries of the density matrix, 
Pll! Pll! Pnn^ and pjvAf can be removed as they include the processes with lepton number 
nonconservation in the active sector and helicity-fiip in the sterile sector, suppressed by the 
small Yukawa couplings and small mass to temperature ratio for A^'s [l3[. As a result, the 
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N /' L '-, N 



FIG. 1: The diagrams contributing to the real potential and rate T in the kinetic equation. 

kinetic equation contains the density matrices pnn and Ptvjv where coherent quantum ef- 
fects are essential, and the diagonal parts of pll and pn, describing concentrations of active 
neutrinos and anti- neutrinos. The latter quantities were not considered by ARS. 

Now we are ready to write our kinetic equations. We will do that for L and only, as 
the equations for L and N are given by CP-conjugation. The medium effects follow from 
computation of real and imaginary parts of the two-point Greens functions on Fig. The 
equation for the sterile neutrinos is 

= [H!:At),PNN] - ^{r^(t),P^;v -P^V} + '^TW{t)F\p,,-pZ)FU{t) . (16) 
For the diagonal part of the active neutrino density matrix we have: 

.dpLL 



diag 



^TF[/(t)(p^;v -p?.V)f/nt)Ft _ '-{T^it),p,,-pl\} 



dt 

Here the following notation is introduced: 



(17) 



H!:.tit) = '^U\t)Knf!Kl,Uit), r^(t) = sin0i7^,(t), T^t) = '-^T K.f.Kl , (18) 

where p^^ and p^^^ are the diagonal equilibrium density matrices for sterile and active 
neutrinos respectively, and sin (p ~ 0.02 |l^ is the ratio of the absorptive part of the diagrams 
to their real part. These equations are the basis for the analysis of the time evolution of the 
asymmetries as a function of the parameters of the i/MSM. 



^ Note that the absorptive parts of these diagrams are in fact not suppressed by the Majorana neutrino 
masses, as was stated by ARS, but by the effective mass of the Higgs boson. However, this effect leads to 



the same numerical estimates of the absorption rates as in (see 0)- 
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The equation p6|) with the third term omitted coincides with that of ARS. The first term 
in (fTBj) describes the oscillations of the sterile neutrinos, whereas the second is responsible 
for absorption and creation of them. The first term in (jl7|) describes the transfer of leptonic 
number from sterile neutrinos to the active ones. It comes from the first diagram of Fig. 
and incorporates the change in F due to the presence of the CP-breaking medium, coming 
from non-trivial pnn- The trace of the second term of (fT^ is exactly equal to the trace 
(with a minus sign) of the first term in ()17|) . ensuring the exact conservation of the lepton 
number. The second term in (fT7j) describes the absorptive processes with active neutrinos 
and ensures that the system eventually thermalizes. The third term in (fTT)|l is a counterpart 
of this one and takes into account the transfer of the active lepton number into the sterile 
sector. It comes from the change in F due to the presence of asymmetries in active neutrinos. 

Our equations have a rich CP-violating structure. In addition to the CP-violating phase 
in the matrix, they contain 3 extra phases, 5^ and ai, a2- Note that the Majorana 
phases /?2 do not appear in the limit of the small Majorana masses we are interested 
in. The equations ()l(i|17|l are to be solved with initial condition pll — p'll — 0|i=o (no 
lepton- flavour asymmetry at the beginning) and patat = 0|t=o (no sterile neutrinos in the 
initial state). 

For sufficiently small Yukawa couplings, when all sterile neutrinos are still out of thermal 
equilibrium at the freeze-out of the sphaleron transitions, the solution can be found pertur- 
batively. This is realized provided f^<2x 10^^^ This regime requires the relatively 
light sterile neutrinos, as from the see-saw formula, and with the use of the absolute values 
of the active neutrino masses one gets M2,3 < 12 GeV. For these masses the processes with 
lepton number non-conservation induced by the Majorana neutrino masses are safely out of 
thermal equilibrium. 

For the higher Majorana masses and larger Yukawa couplings the system can be solved 
numerically. It is clear, however, that if both Nc and A''^ are already thermalized at the 
sphaleron freeze-out, no baryon asymmetry is produced, since the production of the Nd is 
highly suppressed and at that time no neutrino can carry the lepton asymmetry. This leads 
to the requirement that all the sterile neutrino masses must be smaller than the electroweak 
scale, exactly in accordance with our definition of the i/MSM. 

Since the aim of our paper is not the analysis of the complete parameter space of the 
i/MSM but rather the existence proof of a possibility of simultaneous explanation of neutrino 
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oscillations, dark matter and baryon asymmetry, we will consider here the perturbative 
regime only (we checked, however, that numerical solution coincides with the perturbative 
one in the limit of its validity). 

In the region of parameter space where the dark matter constraint is satisfied the gen- 
eration of the lepton asymmetry occurs as follows. First, the sterile neutrinos are created 
in a CP-invariant state (the Jarlskog determinant is very small for K^l) due to the second 
term in eq. (fT^ . However, their density matrix is non-trivial and contains non-diagonal 
terms. Because of that and because of the presence of CP- violating phases in K^, this leads 
to generation of CP-asymmetries in the active neutrino flavours, due to the first term in 
(|17|). Finally, this asymmetry is partially transferred into the total active (or sterile, with 
the opposite sign) asymmetry. 

Now we are ready to solve eqns. (fTBj) and (fTTj) with the use of perturbation theory in . 
For this end they can be rewritten in the integral form, as 

Pnn = (right-hand-side of eq. (fT6|) ) , (19) 

Jo 

Pll = Pll-^ / dt' fright-hand-side of ea. (|T7|l ) (20) 
Jo 

and then solved iteratively. We will present below the leading contributions only and refer 
for the details of calculation to [2^ . 

Let us first discuss the asymmetries in active neutrino flavours induced by the non-thermal 
density matrix of sterile neutrinos. This effect comes from the first term in ()17|) and results 
in a second order effect in A^, = p^.^, — pn: 

^ = ^ T diag (Fp^^Ft - F*p^^F^) . (21) 

To leading order in perturbation theory the elements of pnn come from eq. (pi)|) by neglecting 
the third term on the right-hand side. 

As in the ARS analysis, the asymmetries are produced most effectively at the typical 
temperature Tl ~ (AM^Mq)^/^ and take constant values in the later evolution. Then, when 
Tl ^ Tw ~ 100 GeV, we obtain the following expression 

I _ 7rtsin^0 ^ MqI 

where the effective CP violation parameter is given by 

6fj = Im {F^F)jjFl] . (23) 
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It can be seen that the total lepton number of active neutrinos is zero, i.e. ^fj + ^/j + f^Jj = 0. 
Since the mixing angles sru and sris are very small, as shown in Q, we can safely put 
them to zero in what follows, sru = sri^ = 0. In this case we obtain 

rA 1 I ^ /2/3(/|-/|)Mol 

48-331(1) (AM|2)3 

where a" are given by 

a"" = -4:Sr23Cr23 [sluSlisClis sm{5L + 02)] , 

= 4:Sr23Cr23 [SL12SL13CL13SI23 sm{6L + a2) - cli2Cli3Sl23Cl23 siu ^2] , 

= 4:Sr23Cr23 [SL12SL13CL13CI23 sm{SL + a2) + cli2Cli3Sl23Cl23 siu ^2] • (25) 

It is clearly seen that the phases 6^ and 02 in the mixing matrix are crucial for the 
asymmetry in the active neutrino flavours^. Furthermore, it should be noted that the ob- 
served large mixing angles in the atmospheric and solar neutrino oscillations indicate that 
these parameters are not suppressed but can be 0{1). As an illustration of numerics, we 
shall take the typical values of the Dirac Yukawa couplings to be 



2 J^solMa 2 "^atmMa 
^2 - ^2 ' /3 - -2 ' ^^^) 



where msoi = a/Attt^ ~ 9.1 meV and matm = a/ Am^^m — 51 meV. Then the asymmetries 
are: 

_ a 7r2sin^0 ml^^ml^^Mj Mpi 

\L^L)aoi\Tw — ^ -1-/5^ . , . , ,0 X 2 



48 ■ 33r(|) (AM| 

This equation shows that sizable flavour asymmetries can be generated when the heavier 
"bright" and "clear" neutrinos are highly degenerate in mass. In fact, these flavour asym- 
metries are partially converted into the baryon asymmetry of the universe even if the total 
epton asymmetry is zero, due to effects associated with the masses of the charged leptons 
Iq (for an exact computation of the mass corrections see j3]). However, the conversion 



^ The CP-violating phases which can be found experimentaUy from active neutrino osciUations depend also 
on P2 and thus cannot be used to fix uniquely Sl and a2- 
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2 

rate is suppressed by the factor ~ and asymmetries of (A2,)rr ~ C(O.Ol) are required 

to account for the observed baryon asymmetry, if this mechanism is used. 

However, this is not the end of the story. The above asymmetries in the active neutrino 
flavours trigger the generation of asymmetries in the sterile neutrino sector in a very efficient 
way. This phenomenon is due to the third term in Indeed, the asymmetries An = 

Pnn — Pnn are generated as 



d A AT sin ( 



■T{F^AlF). (28) 



dt 

Since takes constant values after the production time, this equation can be easily solved. 
At the electroweak temperature the asymmetries in the sterile neutrinos are 

(A.)//|T. = ^^(FtA,|,,F),,, (29) 

o J-W 

where A^^Itj^ is defined in (j221)- We should stress that these asymmetries are generated at 
higher order in compared with Al (third versus second order); however the production 
processes continue below the temperature and thus receive a huge enhancement factor ^ 
of the order of Mq/T]^. Since the "dark" sterile neutrino only possesses small Yukawa 
couplings as shown in ((7j), the asymmetries are generated in (A7v)22 and (AAr)33. 

Let us estimate the trace of An, which is crucial for the baryon asymmetry of the universe, 
in the specific parameter choice discussed above. We find from ()27|) that 

I 7rfsin=^0 |2 mi.mi^Mj Mj MqI 

^ ^^T^ .^T.(AM|,)I • ^''^ 

Further simplifications come about as we notice that the successful asymmetry generation 
requires a significant mass degeneracy in M2 and M3 (see below). In order to explain 
the observed active neutrino mass pattern in the normal hierarchy, considered below for 
numerical estimates, one should have a certain hierarchy in Yukawa couplings ^ 
(^ /i^), which gives at the leading order 



J2a''\F^,f = f3'6cP, (31) 



^ Note that when Mq/Tw 00 the perturbative approximation breaks down and one must solve the 
equations exactly. It is clear, however, what happens in this limit: the system will eventually equilibrate 
because of the second terms in ()16I17|I and all asymmetries will go away. However, before this time the 
effect is accumulating, as all sterile neutrinos are out of thermal equilibrium. 
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where 



Sep — 4s_R23C_R23 



Sl12Sli3Cli3((423 + 423)413 " ^13) " ^H^L + 

+ CLi2cii3SL23CL23 (^23 " ^123) " ^2 ■ (32) 



We can see that the CP violation parameter 5c p can be 0{1). Finally, the total asymmetry 
in the sterile neutrino sector is given by 



1 _ 

7r3 sin^(/) m^^im^trnMo Mcf Mn5 



384-351(1) i;6TvK(AM|2)i ^ ^ 

The total lepton number conservation tells us that TrAArlr^^ + Ti A L\Ty^, = and, 
therefore, that the lepton asymmetry in the active neutrino sector is generated ^. It is 
partially converted into the baryon asymmetry due to the rapid sphaleron conversion as 
AB = -f TrA^I^^^ = +|TrA^|^^ Q. 



17|. This completes the computation of the 



baryon asymmetry of the universe in the z/MSM satisfying the dark matter constraint. Ac- 
counting for the entropy factor at the electroweak temperature, the baryon to entropy ratio 
is obtained as 

^...0-^...,.„...0-.»..(^)'(^)\ (34) 

This shows that the correct baryon asymmetry of the universe ^ ^ (8.8 — 9.8) x 10~^^ is 
generated when the heavier sterile neutrinos with the masses, say, 10 GeV are degenerate to 
one part in 10^. This looks like a strong fine tuning but may also indicate the intriguing mass 
relation \M^ — M2I ~ Mi. In this mass range the decays of heavier sterile neutrinos induce 
no significant entropy dilution. The leptogenesis temperature for degenerate neutrinos is 
rather low, Tl ~ 10^ GeV. Note that with this low temperature the right-handed electron is 



16 



3 



can be used. Moreover, 



always thermalized |26| so that the equilibrium formulas of 
if the inverted hierarchy pattern is chosen for the active neutrino masses, even stronger 
degeneracy is needed, since an extra suppression of the baryon asymmetry is coming from 
the small difference between /2 and f^, required in this case. 



^ Equally, the generation of a net lepton number in the sector of active neutrinos can be understood as 
follows. The asymmetries (122(1 in active neutrino flavours are the subject of dissipation described by the 
second term in (|17() . Since does depend on neutrino flavour, these asymmetries evolve differently, 
leading to the total lepton asymmetry. 
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Conclusions. — Let us summarize the obtained results. The i/MSM with three right- 
handed neutrinos with masses smaller than the electroweak scale is the simplest and the most 
economical extension of the Minimal Standard Model. It shares with the MSM its advantages 
(renormalisability and agreement with most particle physics experiments) and its fine-tuning 
problems (the gauge hierarchy problem, fiavour problem, etc). However, unlike the MSM, 
the z/MSM can explain simultaneously three different phenomena, observed experimentally, 
namely neutrino oscillations, dark matter, and baryon asymmetry of the universe. The 
parameter-space of the model is rather constrained (the heavier neutrinos are required to 
be quite degenerate in mass, the Yukawa coupling /i and mixing angles Oru and 9ri3 must 
be very small), which, however, makes this model experimentally testable. The analysis of 
possible experimental signatures of the i/MSM goes beyond the scope of the present letter^. 
We would just like to mention that on the astrophysical side, the best signal would be the 7 
radiation coming from the decay of dark neutrinos. On the particle physics side, as all the 
sterile neutrinos are relatively light, one could imagine that all parameters of the z/MSM, in 
particular the CP-violating phases, will be determined one day (this is clearly a very hard, 
if not impossible experimental challenge because of the smallness of the Yukawa couplings). 
One would then be able to test the baryogenesis formula and in particular its sign. 
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